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On ■ Abstract. Quantum storage and retrieval of light in ion Coulomb crystals using 

cavity electromagnetically induced transparency is investigated theoretically. It is 
' found that, when both the control and probe fields are coupled to the same cavity 

mode, their transverse mode profile affects the quantum memory efficiency in a non- 
trivial way. Under such conditions the control field parameters and crystal dimensions 

■ that maximize the memory efficiency are calculated. 
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^! 1. Introduction 

(N 

Motivated by applications in the field of quantum information processing [H |3l H] , 
quantum memory devices are being investigated in a variety of physical systems and 
^ ■ with different techniques and protocols (for a review see e.g. O [6]). Among the various 

■ criteria used to evaluate the performance of a quantum memory, are generally of interest 

its fidelity, its efficiency, its storage time and its multimode capacity [5]. For optical 
quantum memories, in which an input light pulse is stored into a material system and 
subsequently retrieved, the efficiency can simply be defined as the ratio between the 
energies of the output and input pulses. For an important class of optical quantum 
memories based on electromagnetically induced transparency (EIT) processes in atomic 
ensembles [3 [HI HI [101 [HI [I2l [ISl [IH [15] , the efficiency crucially depends on the optical 
depth of the ensemble [5]. Enclosing the atomic medium in an optical cavity allows for 
substantially increasing the effective optical depth experienced by the light and cavity 
EIT protocols have been proposed to achieve high storage efficiencies [Ml [TTj [TU |19l [20] . 

Among the variety of atomic systems studied in connection with quantum memories, 
an ion Coulomb crystals positioned in an optical cavity has been suggested as a good 
candidate to realize a high-performance quantum memory, potentially meeting the 
criteria mentioned above [21]. Recently, this analysis has been backed up by a series of 
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key experimental results. Substantial effective optical depths can indeed be realized by 
strongly coupling large ion Coulomb crystals to a cavity field and long coherence times 
can be achieved in such a system |22l |23]. Furthermore, strong coupling to various 
spatial cavity modes has been demonstrated [21], which is promising for multimode 
storage. Finally, cavity EIT has recently been observed with ion Coulomb crystals 
using an a//-cavity geometry in which both the probe and control fields are coupled to 
the same cavity mode [25] . 

In this specific cavity EIT geometry, it has been observed, both theoretically [26] 
and experimentally [25], that the spatial transverse profile of the control field has a 
significant effect on the probe transmission/reflection lineshapes and dynamics. This 
is in contrast with more standard cavity EIT configurations in which the control 
field is free-propagating and has typically a much larger extension than the probe 
field [271 EHl EHl [301 SB ESI EH]. Since the all-cavity geometry is one natural realization 
of the scheme presented in [21], it is interesting to investigate its implications for 
the performance of such an ion Coulomb crystal-based quantum memory. In this 
paper, we extend the existing theoretical models for cavity EIT-based light storage and 
retrieval [161 113 113 UHl EO] to the all-cavity configuration and numerically investigate 
the effect of the transverse mode profile of the fields on the quantum memory efficiency. 
We find in particular that the optimal efficiency depends not only on the cooperativity 
parameter, but also on the radial extension of the crystal, as a result of the more complex 
spatial mode structure defined by the fields inside the atomic medium. Using parameters 
taken from current experiments with ion Coulomb crystals [221 ES], our simulations 
predict that similarly high-efficiencies (> 90%) should however be obtainable in the 
all-cavity configuration. 

The paper is structured as follows: in sec. [2] the theoretical model for the light 
storage and retrieval is presented, starting with a general description in sees. 12.11 and 
12.21 an optimization of the quantum memory efficiency in the standard configuration 
for smoothly- varying single-photon input pulses, following the approach of ref. [20] , 
in sec. 12.31 and the "cylindrical shell" model used for the simulations of the all-cavity 
configuration in sec. 12.41 Section |3] presents the results of the numerical simulations 
based on typical experimental parameters for ion Coulomb crystals (sec. 13. ip and for 
light fields coupled to the cavity fundamental TEMqo mode (sec. 13. 2p and a higher-order 
Laguerre- Gauss LGqi mode (sec. 13. 3p . A brief conclusion on the implications of these 
results for experiments are given in sec. ID 

2. Theoretical model 

2.1. Description of the model 

Inspired by the experiments of [25] we base our description of the light-matter interaction 
on the model developed in [26]. We consider an ensemble of three-level A atoms with 
two ground states, |1) and |2), and an excited state |3). The atoms interact with a cavity 
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probe field on the |1) — )■ |3) transition and a classical control field Q on the |2) — > |3) 
transition. Both fields are assumed to be resonant with the atomic transitions and 
the cavity is tuned to resonance with the probe field. Since we are mostly interested 
in discussing the effects of the transverse mode profile of the fields, we neglect the 
longitudinal variation of the atom-field couplings. This goes for atoms in a running- 
wave cavity [291 |32l [31] , but can also apply to a standing- wave cavity geometry for atoms 
with either well-defined positions with respect to the cavity standing- wave [281 EH 126] or 
atoms whose motion average out the longitudinal standing-wave structure of the fields 
over the relevant timescales [251 126]. The atomic ensemble extension is also assumed to 
be much smaller than the cavity field Rayleigh range. The cavity is furthermore taken 
to be single-ended and with lossless mirrors (fig. [1]). 




Figure 1. (Color online) Left: Three-level A atomic structure considered. Right: a 
cylindrically symmetric atomic medium, composed of A atoms, is enclosed in a single- 
ended optical cavity where it interacts with a probe field and a control field in an EIT 
situation. Two configurations are investigated: one in which both field are coupled 
to the same cavity mode - and have therefore the same waist, and one in which the 
control field has a constant intensity profile over the ensemble. 

We will also assume that the fields injected into the cavity are smooth pulses whose 
envelopes are slowly varying with respect to the cavity field and atomic decay rates, in 
the sense of [20] . For the light storage and retrieval dynamics we will follow the approach 
of refs. [m iini |20] and assume a standard write-store-read temporal sequence. During 
the write phase the input probe field is injected into the cavity and the control field 
is adiabatically turned off to ensure temporal impedance matching [18]. Both fields are 
turned off during the store phase, and the control field is adiabatically turned back on 
during the readout, causing the emission of an output probe pulse. 

Denoting by ajl'J = |/i)(z/|fc the individual atomic operator for the k—th atom 
positioned at and by a the annihilation operator for the intracavity probe field, the 
interaction Hamiltonian in the rotating wave approximation and the rotating frame 
reads 

H = -hgJ2 '^p{rk)4i(^ - ^^(t) J2 '^c{rk)42 + h-c (1) 

k k 

where g is the single-atom maximal coupling rate (at the center of the cavity mode). 
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'^p,c{f) are the probe and control field transverse mode profiles (with the longitudinal 
dependence neglected) and Q{t) the time- varying Rabi frequency of the control field. 

2.2. Equations of motion in EIT 

We consider a typical EIT regime in which all the atoms are initially in level |1) and 
in which the control field is much more intense than the probe field, assumed to be 
at the single-photon level. One can then assume that almost all the atoms stay in |1) 
at all times and perform a standard first-order treatment in the probe field [2]. The 
Heisenberg equations of motion for the relevant operators, namely the intracavity probe 
field a, the atomic optical coherences a^^^ and the atomic ground state coherences i 
are given by 

a = -Kfl + ig^ '^p{rk)dr[^^ + V2Kain, (2) 

k 

^13 = -l^is + W'^p{rk)a + in{t)^,{r,)a['^ + F^K (3) 
^12 = -70^1? + tn*m,ir,)a[l^ + fS\ (4) 
and the input-output relation 

Clout = V2K,a - Clin, (5) 

where k is the cavity field decay rate, and dout are the annihilation operators 
associated with the input and output probe fields, respectively. 7 and 70 are the atomic 
dipole and ground state coherence decay rates, respectively, and F^^^ and F^2^ the 
corresponding Langevin noise operators. 

Following ref. [20], we assume a single-photon input and calculate the quantum 
memory efficiency by first solving in time the semiclassical counterparts of ©-([S]) for 
given input probe and control field pulses 

a = -na + ig^ -^p{rk)crfj + v^a^^, (6) 



!? = -l^u + ^9'^p{rk)a + zfi(t)^e(rfe)aiy , (7) 
^i? = -70^!? + tn\t)^>,{n)a'S, (8) 
o-out = \/2Ka — a^n, (9) 
and, secondly, by computing 

_ j,\aout{t)\^dt 

where the subscripts w and r refer to a summation over the whole duration of the write 
and read process, respectively, rjtot then represents the ratio of the number of retrieved 
photons to the number of incoming photons, and provides a good measure of the quality 
of the mapping (for other measures see e.g. O |6]). To simplify the discussion we will in 
the following neglect the decay of the ground state coherence during the whole process 
and set 70 = 0. The duration of the storage phase is then simply chosen such that the 
dynamics of the write and read phases occur in well-separated time windows. 
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2.3. Efficiency optimization in the case of an extended control field 

If the waist of the control field is much larger than the probe field, as in many EIT 
experiments, one can neglect the transverse variations of the control field intensity 
over the section of the atomic ensemble that interacts with the probe field. Following, 
e.g. [26], one can define an effective number of atoms interacting with the probe field 

N = Y,^,{nf (11) 

k 

and collective operators for the ground state coherence and the optical dipole by 



S 



Equations ([6])- ([9]) can be straightforwardly rewritten as 

a = —Ka + igNP + V^Kcun, (13) 
p = + ig^a + iVL{t)S, (14) 
S = tQ*{t)P, (15) 

Oout = V2Ka - Oin, (16) 

where = gVN is the collectively enhanced coupling rate [2]. 

In this collective mode picture, the input photonic state is mapped during the write 
phase onto a collective spin-wave described by S and one can define a write efficiency by 
taking the ratio of the number of atomic excitations and the number of input photons, 

|ai„(t)|2rft' 



where is the end time of the write phase [20]. Similarly, one can define a readout 
efficiency by the ratio of the number of output photons and the number of atomic 
excitations before readout 

where is the start time of the read phase. In the adiabatic limit, i.e. for input pulses 
with duration T such that 2TC7 ^ 1, one can derive the optimal control pulse that 
maximizes the read and write efficiencies, which can be shown to scale as [19l |20] 

2C 

where 



1 + 2C" ^^^^ 



C^(^ (20) 

is the cooperativity parameter |j. In absence of decoherence during the storage phase, 
the optimal total efficiency thus scales as 

(21) 

I Note the factor 2 difference with respect to j20| . 
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and increases with the effective number of atoms defined by the spatial overlap of the 
ensemble and the probe field [Eq. f|TT]) ]. 

2.4- Effect of the control field's transverse profile 

If the waist of the control field Wc is no longer very large, but comparable to that of 
the probe field Wp - as it will be the case in an all-cavity geometry for instance [25| [26] 
- the previous results no longer apply, and one must evaluate the effect of the control 
field transverse profile on the storage and retrieval efficiencies. With the approximations 
made in Sec. 12.21 and having specifically in mind ion Coulomb crystals as the physical 
storage medium, we assimilate the atomic ensemble to a cylinder with length L and 
radius R (fig. [1]), and slice it into n cylindrical shells of thickness d <^ Wp, Wc {R = nd). 
We also assume that d is much larger than the mean interparticle distance. Although it 
is not essential, we take the atomic density p to be constant throughout the ensemble, 
which is the case for large ion Coulomb crystals in linear Paul traps , and we consider 
cavity modes with cylindrical symmetry. We proceed by defining collective operators 
for the jth-slice as 

Pj = aUjCTi^ ^ Sj = anja^2 } (22) 

where the subscript j refers to an atom in the j-th slice with position rj = d{j — 1/2) 
{j = l..n) and where a = pL is the atomic cross-sectional density. The corresponding 
semiclassical equations of motion are then 

n 

a = -na + ig ^ '^p{rj)Pj + V2K,ain, (23) 

Pj = -^Pj + i gauj -^p {rj )a + in{t)-^c {rj ) 5"^ , (24) 
S,=tn\t)<if,ir,)S,, (25) 

where the modefunctions \l/p,c are evaluated at rj. 

It is clear from eqs. that, unless |\l/c('^j)| = 1 like in the previous 

section, it is no longer possible to define collective spatial eigenmodes of the problem that 
would yield closed equations of the form f[T^ . f[T^ . f[T^ . In particular, these equations 
show that adjacent shells are coupled together by the control field. The spatial mapping 
of the probe field onto the ground state spin now depends on the control field transverse 
profile, in addition to that of the probe field. Because of the intershell coupling during 
the mapping the radial extension of the ensemble now becomes a parameter which affects 
the memory efficiency in a non-trivial way. 

In the following section we numerically solve these equations of motion for a fixed 
probe field pulse basing ourselves on the analytical control field pulse derived from 
the temporal optimization of sec. 12.31 in absence of effects due to the control field 
transverse profile. Note that, because of the impossibility of defining analytically a 
spatial collective spin-mode during the write or read phase, and thereby of defining 
write or read efficiencies, this optimization is performed numerically using the total 
efficiency rjtot as a figure of merit. 
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3. Numerical results 

3.1. Physical system considered and input parameters 

To solve the problem of optimizing the quantum memory efficiency in the conditions of 
the previous section we take for the physical storage medium an ion Coulomb crystal, 
trapped and laser-cooled in a linear Paul trap, with the optical cavity positioned 
along the trap axis, as in [22l [241 EHl [25]. Although single-component ion Coulomb 
crystals (i.e. consisting of only ion species) have spheroidal shape and may therefore 
deviate from the cylindrical shell model (unless they are sufficiently prolate), the inner 
component of a prolate two-species crystal can be assimilated to a good approximation 
to a uniform density cylinder [371 [M]- We consider the cavity EIT configuration 
used in [25], in which both the control and probe fields are frequency degenerate 
and orthogonally polarized in order to create EIT between Zeeman sub-states of the 
3d ^D^i2 sub-level of ^"^Ca"*". In these experiments both fields are coupled resonantly or 
near-resonantly to the same cavity mode (TEMqo). We thus take Gaussian transverse 
profiles "^p^cij) = 6xp(— r^/w^^), and compare the extended control field configuration 
{wc — )■ oo) and the finite control field waist configuration {wc = Wp). 

For the 11.8 mm-long, close to confocal cavity of [25] with an incoupling mirror 
transmission of 1500 ppm and an interaction on the 3d ^0^/2, fnj = +3/2 — )■ 
Ap ^Pi/2,mj = +1/2 (probe) and 3d ^Dy2,mj = -1/2 ^ Ap ^Pi/2,mj = +1/2 
(control) transitions, one finds {g, k, 7) = 27r x (0.37, 1.5, 11.3) MHz[§|. With a radius of 
curvature of 10 mm, the waist of the probe field at the center of the cavity is Wp = 37 
/xm. For crystals with typical length of a few mm and radius of up to a few hundreds of 
microns, neglecting the longitudinal curvature of the fields over the crystal modevolume 
is well-justified [221 l2l]. 

For the sake of the discussion, we assume for the probe field a hyperbolic secant 
input pulse of the form 



where T is the probe pulse duration [jjj. In the adiabatic limit (TC7 ^ 1) considered 
previously in sec. 12.31 and the extended control field configuration {wc — >■ 00), the control 
field pulses which optimize the write and read efficiencies are given by [20] 



account for the longitudinal averaging over the standing-wave structure. 

II This form for the probe pulse is taken for convenience, as one gets an analytical expression for the 
optimal control field pulse [ISi, but the numerical simulations show that the exact form of the input 
pulse is not critical. 



ain{t) = — sech(2t/T), 



(26) 




(27) 



for the write phase, and its time-reversed counterpart 




8 



for the readout phase. In the extended configuration the pref actor A is equal to unity. As 
we will see in the next section, the previous control field temporal profiles are still found 
to be optimal with respect to maximizing the efficiency in the finite waist configuration, 
the main difference being in the optimal control field amplitude scaling factor A. 

3.2. Results for the TEMqq mode 



-T- 




-'^O 10 20 30 

tins] 

Figure 2. Temporal storage and retrieval sequence for (a) an extended control field 
(wc oo), (b) a control field with finite waist {wc = Wp) and non-optimized amplitude 
[A = \) and (c) a control field with finite waist {wc = Wp) and optimized amplitude 
{A = 2.45). The dotted and solid curves are the normalized input and output field 
intensities, |aj„(i)p and \aout{t)\'^, respectively. The dashed curve shows the control 
field intensity ^(t)^, scaled such that its maximum value is 0.5. = = 10 fis. See 
text for the value of the other parameters. 

Figure [2^ shows the results of a storage and retrieval sequence for a crystal with 
density p = 6.1 x 10^ cm~^, length L = 3 mm and radius R = 100 fim, an input probe 
pulse of the form fl2B]) with duration T = 2 fis and an extended control field of the 
form (127|) . (128!) with A = 1. For such a large crystal {R ^ Wc), the effective number 
of ions as defined by (fTTl) is = 3936, yielding a cooperativity parameter C ~ 16.7. 
The write and read efficiencies are found to be rj^ = 0.970 and rjr = 0.971, respectively, 
close to the theoretical value of r]"^^ = 0.971, and yielding a total efficiency rjtot = 0.942 
iVtot = 0.943). Figure [2b shows the results of the same sequence and parameters for a 
control field with Wc = Wp. It is clear that, during the write phase, perfect temporal 
impedance matching is not achieved, as a substantial amount of the incoming light is 
refiected. This is in itself not surprising as the prefactor A = 1 is only optimal in 
the extended configuration, and one expects that, in the finite waist configuration, the 
ions see on average a control field with lower Rabi frequency. As can be seen from 
fig. [2]3, close to perfect impedance matching can be recovered by increasing the control 
field amplitude {A ~ 2.45). The total efficiency rjtot = 0.667 remains, however, lower 
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than in the corresponding extended configuration. We checked numerically that varying 
the control field pulse switching time and shape, or having different amplitudes/time 
evolutions profiles during write and read, does not increase the efficiency. We also 
checked that these results do not significantly depend on T as long as one stays in the 
adiabatic limit. 




0.5 1 1.5 2 2.5 

R/ w 



P 

Figure 3. Variation of the total efficiency rjtot as a function of the crystal radius R, for 
the same parameters as in fig. [2] The dashed line shows the extended control field case 
while the solid line shows the finite waist case for which the control field amplitude has 
been optimized for each radius. The inset shows the variation of the effective number 
of ions [eq. ([TT])] as a function of the crystal radius R. 

The lack of critical dependence of the optimal efficiency with respect to the temporal 
parameters found in the simulations seems to indicate, in agreement with the analysis 
of sec. 12.41 that the spatial profile of the control field now plays a significant role in the 
mapping process. To investigate this effect further we show in fig. [3] the variation of 
the effective number of ions and the (temporally) optimized efficiency as a function of 
the crystal radius in the two configurations. The other parameters are kept the same 
as previously. While in the extended control field configuration the efficiency increases, 
together with A^, with the crystal radius and saturates when R ^ Wp, it reaches a 
maximum for R Wp in the finite waist configuration, before decreasing and reaching a 
constant level at high radius. The total efficiency peaks at ~ 0.95wp with a value of 
0.914 {A = 1.5), closer to the theoretical value of 7]°^^ = 0.932 (for this radius = 3279 
and C = 13.9). The decrease for R > Wp radius may appear surprising, since more ions 
are being added to the crystal and one could expect an enhanced efficiency due to the 
stronger coupling to the probe field. However, for a finite waist control field, the spatial 
spin mode defined during writing - or, for that matter, during reading - is no longer 
that of the probe, but depends on the overlap of both the probe and the control field 
transverse profiles in the crystal. As the radius of the crystal is increased the stored 
photonic excitation is spread more and more over shells with higher radius. The spatial 
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atomic mode then no longer resembles the spatial mode of the probe, which leads to 
a decrease in efficiency in the writing process. A similar phenomenon then takes place 
in the reading process, as the atomic spatial excitation proffie is no longer optimally 
matched to the spatial light mode profile. 
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Figure 4. Variation of the radial density of excitations after writing in the the 
extended (dashed) and finite waist (sohd) configurations, for different crystal radii 
[(a) R = 0.5wp (b) R = O.OSuij, (c) R = 2.1w^. Other parameters as in Fig. [3] 



To illustrate this phenomenon, we show in fig. H] the variation of the radial density 
of excitations after writing, which is proportional to the surface probability of finding 
the photonic excitation in the j-th shell, at the end of the write phase, again in both 
configurations and for crystals with different radii. The radial density of excitations 
s(rj) is obtained by normalizing the squared modulus of the j-th shell operator mean 
value Sj{Ty^) by the number of ions in the shell rij, 

sir,) = \S,iT^)\/n,\' (29) 

In the extended configuration, one sees that s{rj) reproduces well the spatial Gaussian 
mode profile of the probe field, as expected from the analytical predictions of sec. 12.31 In 
the finite waist configuration, the spatial atomic mode defined by s{rj) has reasonable 
overlap with the ideal extended configuration mode (i.e. the probe field mode) for 
R < Wp, but clearly deviates from it as the crystal radius increases and the coupling 
between the shells causes the excitation to spread more and more radially into the 
crystal. 

To summarize the results and show in particular that the decrease in efficiency 
for large radii is always substantial in the finite waist configuration, regardless of the 
length or density of the crystal, fig. [5] shows the variation of the optimized efficiency as a 
function of the crystal dimensions, L and R, for a density of 6.1 x 10^ cm~^. The range 
chosen for the dimensions is typical of current experiments with ion Coulomb crystals 
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Figure 5. (Color online) (a) Variation of the cooperativity as a function of the 
crystal dimensions L and R. (b) Variation of the optimized efficiency in the extended 
configuration versus L and R. (c) Variation of the optimized efficiency in the finite 
control field waist configuration (wc = Wp) versus L and R. The crystal density is 
p = 6.1 X 10® cni^'^ and the pulse duration T — 2 /is. 



in cavity [22l [38] . The optimized quantum memory efficiency in the extended control 
field configuration is found to agree well with the predictions from the analytical model 
[eq. [21], and is an increasing function of the crystal length for all radii. This is no longer 
true in the finite waist configuration for which an optimal radius - of the order of Wp - 
exists, for all lengths. However, it can be seen that by choosing the radius of the crystal 
appropriately to "match" the waist of the cavity mode one can achieve similarly high 
quantum memory efficiencies (> 90%) as in the extended configuration. 

3.3. Higher-order modes 

In this last section we investigate these spatial mode effects on the storage using higher- 
order spatial cavity modes. On the one hand, this is motivated by the fact that 
collective strong coupling with higher-order cavity modes has been demonstrated using 
ion Coulomb crystals [23], which is promising for multimode (spatial) storage. On the 
other hand, in view of the previous results, one can wonder how the conclusions drawn 
for the TEMoo mode generally hold for higher-order cavity modes, and in particular, if 
some modes are less sensitive to these spatial effects. A general analysis is beyond the 
scope of the present paper and we will only focus in this last section on the case of a 
first-order Laguerre- Gauss cavity mode LGqi, which preserves the cylindrical symmetry 
of the problem. We thus assume that the probe field radial modefunction is now given 
by \E'p(r) = \/2{r/wp)exp{—r^/Wp). We then compare the quantum memory efficiency 
in an extended control field configuration and in a configuration where the control field 
has the same transverse profile as the probe field (\E'c(^) = ^pl'"))- 

Figure [6^ shows the variation of the effective number of ions defined by eq. ([TTj) as 
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Figure 6. (Color online) Variation of the quantum memory efficiency rjopt as a function 
of the crystal radius R, for the TEMqo (black) and LGqi (green) modes and for a 
crystal with L — 3 mm and p = 6 x 10^ cm~^. The dashed lines show the extended 
configuration and the solid lines the finite waist configuration. The inset shows the 
corresponding variations of the effective number of ions [eq. pTjl ] as a function of R, 
for both modes. 

a function of the crystal radius, for a crystal with fixed length and density (L = 3 mm 
and p = 6.1 X 10^ cm^^). As expected, N increases less rapidly at small radii than for 
the fundamental mode, on account of the lower coupling at the center of the mode, but 
saturates at the same value for large radii, because the orthonormal character of the 
mo defunct ions. In fig. [Ub are represented the corresponding variations of the quantum 
memory efficiency, both for the TEMqo and LGqi modes and in the extended and finite 
waist configurations. If a similar behavior is qualitatively observed for the LGqi mode, 
one can see that the optimal radius is larger for the LGqi mode {R ^ 1.35wp) than for 
the TEMqo mode [R ~ 0.95wp). This can be explained by the fact that the regions of 
high radial intensity for the fields are now located farther away from r = 0, so that the 
spreading of the photonic excitations into large radius shells occurs at larger R and is 
less pronounced. One can carry out a similar analysis as for the fundamental mode and 
calculate the variation of the radial density of excitations after writing as a function of 
the crystal radius. The results are represented in fig. [7] and show the same qualitative 
conclusions as drawn previously in the case of the fundamental mode. 

4. Conclusion 

We have investigated the efficiency of a cavity EIT-based quantum memory in which 
both control and probe fields are coupled to the same cavity mode. Due to the complex 
spatial atomic mode defined during the EIT process between both fields during the 
write and read phase of the memory, the optimal efficiency is found to depend not only 
on the cooperativity parameter, but also the crystal radius. Using parameters from 
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Figure 7. (Color online) Variation of the mean surface number of excitations after 
writing in the the extended (dashed) and finite control field waist (solid) configurations, 
for different crystal radii [(a) R — O.bWp (b) R — l.SSwp (c) R = 2.7wp] and for the 
LGoi mode. Other parameters as in Fig. [S] 



current experiments with ion Coulomb crystals [221 EH], our simulations predict that 
high-efficiencies (> 90%) should however be obtainable in this specific configuration. 

This theoretical investigation clearly implies that an experimental realization of 
such an all-cavity EIT quantum memory based on a cylindrical ion Coulomb crystal 
can be optimized by choosing a crystal with radius matching the waist of the cavity 
modes. For single species crystals as used in |22l [25] , the optimum crystal radius has 
been found to differ slightly because of the spheroidal (non-cylindrical) shape of such 
crystals [37], but similar conclusions hold. For two-species crystals, the lighter species 
typically takes the shape of a nearly perfect cylindrical rod surrounded by the other 
ion species [371 EH |39], and matches perfectly the situation considered theoretically. 
Applying isotope selective photoionization [lO], Coulomb crystals consisting of two 
calcium isotopes can easily be created with varying compositions and shapes [381 EH] > 
providing an ideal situation to test the predictions and optimize the storage conditions. 
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